Hubbard Fermi surface in the doped paramagnetic insulator 
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We study the electronic structure of the doped paramagnetic insulator by finite temperature 
Quantum Monte-Carlo simulations for the 2D Hubbard model. Throughout we use the moderately 
high temperature T = 0.33t, where the spin correlation length has dropped to < 1.5 lattice spacings, 
and study the evolution of the band structure with hole doping. The effect of doping can be best 
described as a rigid shift of the chemical potential into the lower Hubbard band, accompanied by 
a transfer of spectral weight. For hole dopings < 20% the Luttinger theorem is violated, and the 
Fermi surface volume, inferred from the Fermi level crossings of the 'quasiparticle band', shows a 
similar doping dependence as predicted by the Hubbard I and related approximations. 
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Since the pioneering works of Hubbard |il| , the metal- 
insulator transition in a paramagnetic metal has been 
the subject of intense study. Despite this, our theoret- 
ical understanding of this phenomenon is quite limited. 
Hubbard's original solutions to the problem, the so-called 
Hubbard I-HI approximations, have recently faced some 
criticism. One fact which is frequently held against his 
approximations or the closely related two-pole approxi- 
mations 0-0] is the difficulty to reconcile them with the 
Luttinger theorem. This can hardly be a surprise in that 
all of these approaches rely on splitting the electron cre- 
ation operator into two 'particles' which are exact eigen- 
states of the interaction term Hu: 
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with [HuJl„ ] = Udl^, and [Hc/,^^ ] = 0. The inter- 
action term is therefore treated exactly, approximations 
are made to the kinetic energy. This is precisely the oppo- 
site situation as compared to the perturbation expansion 
in U, which leads to the Luttinger theorem. At half- 
filling the two 'particles' d] ^ and c]^, whose energy of 
formation differs by U, then form the two separate Hub- 
bard bands. The effect of doping in both the Hubbard 
I approximation or the two-pole approximations consists 
in the chemical potential cutting gradually into the top 
of the lower Hubbard band, in much the same fashion 
as in a doped band insulator. On the other hand the 
spectral weight along the lower Hubbard band deviates 
from the free-particle value of 1 per momentum and spin 
so that the Fermi surface volume (obtained from the re- 
quirement that the integrated spectral weight up to the 
Fermi energy be equal to the total number of electrons) 
is not in any 'simple' relationship to the number of elec- 
trons - the Luttinger theorem must be violated. 
In this manuscript we wish to address the question as 
to what really happens if a paramagnetic insulator is 
doped away from half-filling, by a Quantum Monte Carlo 
(QMC) study of the 2D Hubbard model. We use the 



value U/t = 8 and work throughout at the moderately 
high temperature T = 0.33t. This temperature is small 
compared to both the bandwidth, U, and the gap in the 
single particle spectrum (see Figure ||). The main ef- 
fect of T is the destruction of antiferromagnetic order, 
as discussed in our previous paper [||. We therefore be- 
lieve that our study realizes to good approximation the 
situation for which Hubbard's solutions were originally 
designed: a paramagnetic system in the limit of large U, 
at a temperature which is small on the relevant energy 
scales. 

Below, we present results for the single particle spectral 
function and its doping dependence. These data show, 
that the Hubbard I approximation is in fact considerably 
better than commonly believed: the effect of doping in- 
deed consists mainly of a progressive shift of the chemical 
potential fi into the band structure of the insulator. The 
Fermi surface volume, if determined in an 'operational' 
way from the single particle spectral function, indeed is 
not consistent with the Luttinger theorem. 
We start with a brief discussion of the band structure at 
half-filling, see Figure |^, which shows the single particle 
spectral function. We note that this is quite consistent 
with previous QMC work W. For comparison the two 
bands predicted by the Hubbard I approximation. 
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are also shown as the dashed dispersive lines 
{ek——'2t{cos{kx) + cos{ky) ) is the nonintcracting dis- 
persion). These provide at best a rough fit to those 
parts of the spectral function which have high spectral 
weight. Inspection of the numerical spectra shows quite 
a substantial difference between the numerical and the 
Hubbard-type band structures: the latter always give two 
bands, whereas in the numerical spectra one can rather 
unambiguously identify 4 of them, denoted as as i?. A, 
A' and B' (see the Figure). None of these bands shows 
any indication of antiferromagnetic symmetry; together 
with the short spin correlation length [pj 






FIG. 1. Gray-scale plot of the single particle spectral den- 
sity for the 20 x 20 cluster at half-filling. The gray scale gives 
the intensity of spectral weight at the respective k — uj point. 
Also shown are the Hubbard I bands (dashed lines) and the 
4-bands obtained by solving (H) (full lines). 

this shows that we are really in the paramagnetic phase. 
We found that to model this 4-band structure one can 
inroduce two additional dispersionless bands at energies 
of E± = ^ ± e. We now allow mixing between each of 
these dispersionless bands and the respective Hubbard 
band, as would be described by the Hamilton matrix 
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Using the values e — it and V — t, the resulting 4-band 
structure provides at least a qualitatively correct fit to 
the numerical data. We stress that at present we have no 
'theory' for these two additional bands. Equation (0) is 
just a phenomenological ansatz to fit the numerical band 
structure. We note, however, that a 4-band structure 
which has some similarity with our results has recently 
been obtained by Pairault et al. using a strong coupling 
expansion g. 

We do not, however, pursue this issue further but turn 
to our main subject, the effect of hole doping. Figure g 
shows the development of A(fc,a;) with doping. Thereby 
the A{k,uj) for different hole concentrations have been 
overlaid so as to match dominant features, and the chem- 
ical potentials for the different hole concentrations are 
marked by lines. It is quite obvious from this Figure 
that the 2 bands seen at half-filling in the photoemission 
spectrum persist with an essentially unchanged disper- 
sion. The chemical potential gradually cuts deeper and 
deeper into the A band, forming a hole-like Fermi surface 
centered on (tt, tt), the top of the lower Hubbard band. 
The only deviation from a rather simple rigid-band be- 
havior is an additional transfer of spectral weight: the 
part of the A-band near (tt, tt) gains in spectral weight, 
whereas the _B-band looses weight. The loss of the B 
band cannot make up for the increase of the A band, but 
rather there is an additional transfer of weight from the 



upper Hubbard bands, predominantly the A' band. This 
effect is quite well understood Q. The A' band seems 
to be affected strongest by the hole doping and in fact 
the rather clear two-band structure visible near (tt, tt) 
at half- filling rapidly gives way to one broad 'hump' of 
weight. Apart from the spectral weight transfer, however, 
the band structure on the photoemission side is almost 
unaffected by the hole doping - the dispersion of the ri- 
band becomes somewhat wider but does not change ap- 
preciably. In that sense we see at least qualitatively the 
behavior predicted by the Hubbard I approximation. 




FIG. 2. Single particle spectral function for the 8x8 clus- 
ter and different electron density (n) . The chemical potential 
at half filling is the zero of energy, the spectra for different 
(n) are rigidly shifted relative to one another so as to match 
dominant features The chemical potentials for the different 
(n) are given by vertical lines. 

Next, we focus on the Fermi surface volume. Some care 
is necessary here: first, we cannot actually be sure that 
at the high temperature we are using there is still a well- 
defined Fermi surface. Second, the criterion we will be 
using is the crossing of the A band through the chemical 
potential. It has to be kept in mind that this may be 
quite misleading, because band portions with tiny spec- 
tral weight are ignored in this approach (see for example 
Ref. ijlffl for a discussion). When thinking of a Fermi 
surface as the constant energy contour of the chemical 
potential, we have to keep in mind that portions with low 
spectral weight may be overlooked. On the other hand 
the fact that a peak with appreciable weight crosses from 
photoemission to inverse photoemission at a certain mo- 
mentum is independent of whether we call this a 'Fermi 
surface' in the usual sense, and should be reproduced by 
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fore has to be kept in mind that in the following we 
are basically studying a 'spectral weight Fermi surface', 
i.e. the locus in k space where an apparent quasiparticle 
band with high spectral weight crosses the chemical po- 
tential. With these caweais in mind, Figures Hand H show 
the low-energy peak structure of A{k,uj) for all allowed 
momenta of the 8x8 cluster in the irreducible wedge 
of the Brillouin zone, and for different hole concentra- 
tions. In all of these spectra there is a pronounced peak, 
whose position shows a smooth dispersion with momen- 
tum. Around (tt, tt) the peak is clearly above /x, whereas 
in the center of the Brillouin zone it is below. The locus 
in fc-space where the peak crosses fi forms a closed curve 
around (tt, tt) and it is obvious from the Figure that the 
'hole pocket' around (tt, tt) increases very rapidly with 
d. To estimate the Fermi surface volume Vp we assign 
a weight Wk of 1 to momenta k where the peak is be- 
low /i, 0.5 if the peak is right at fi and if the peak is 
above /i. Our assignments of these weights are given in 
Figure 0. The fractional Fermi surface volume then is 
Vp — -^X^k^k, where iV = 64 is the number of mo- 
menta in the 8x8 cluster. Of course, the assignment of 
the Wk involves a certain degree of arbitrariness. It can 
be seen from Figures and 0,however, that our lUk 



FIG. 3. Single particle spectral function for all k-points of 
the 8x8 cluster in the irreducible wedge of the Brillouin zone. 
For each k the weight w^ is given. 
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FIG. 4. Same as Figure B^ for lower electron densities, 
any theory which claims to describe the system. It there- 




FIG. 5. Fermi surface volume as estimated from the single 
particle spectral function, plotted versus the concentration of 
holes in the half-filled band. The dashed line gives the value 
predicted by the Luttinger theorem, W = f ■ 

would in any way tend to underestimate the Fermi sur- 
face volume, so that the obtained Vp data points rather 
have the character of a lower bound to the true Vp ■ Even 
if we take into account some small variations of Vp due to 
different assignments of the weight factors, however, the 
resulting Vp versus S curve never can be made consistent 
with the Luttinger volume, see Figure 0. The deviation 
from the Luttinger volume is quite pronounced at low 
doping. Vp approaches the Luttinger volume for dopings 
~ 20%, but due to our somewhat crude way of determin- 
ing Vp we cannot really decide when precisely the Lut- 
tinger theorem is obeyed. The Hubbard I approximation 
approaches the Luttinger volume for hole concentrations 
of « 50%, i.e. the steepness of the drop of Vp is not 



reproduced quantitatively. The latter is somewhat im- 
proved in the so-called 2-pole approximation |j^. For 
example the Fermi surface given by Beenen and Edwards 
g for (n) = 0.94 obviously is very consistent with the 
spectrum in Figure H for {n) — 0.95. 
In summary, we have studied the doping evolution of 
the single particle spectral function for the paramagnetic 
phase of the 2D Hubbard model, starting out from the 
insulator. As a surprising result, we found that in this 
situation the Hubbard I and related approximations give 
a qualitatively quite correct picture. The main discrep- 
ancy between the Hubbard I and the so-called 2-pole ap- 
proximation and our numerical spectra is the number of 
'bands' of high spectral weight, which is 4 in the numer- 
ical data. This is no reason for concern, because we have 
seen that adding two more bands allows for an quite rea- 
sonable fit to the numerical band structure and one might 
expect that finding a somewhat more intricate decoupling 
scheme for the Hubbard I approximation or a suitable 4- 
pole approximation should not pose a major problem. 
The greatest success of the Hubbard-type approxima- 
tions, however, is a qualitatively quite correct descrip- 
tion of the evolution of the 'Fermi surface'. The effect 
of doping consists of the progressive shift of the chemical 
potential into the topmost band observed at half-filling, 
accompanied by some transfer of spectral weight. The 
Fermi surface volume, determined in an 'operational way' 
from the band crossings, violates the Luttinger theorem 
for low hole concentrations and does not appear to be in 
any simple relationship to the electron density. The Lut- 
tinger sum rule is recovered only for hole concentrations 
around 20%. 

It is interesting to note in this context that a recent 
study of the momentum distribution in the t-J model 
by high-temperature series expansion |ll| ] has also pro- 
vided some evidence for a 'Fermi surface' which encloses 
a larger value than predicted by the Luttinger theorem. 
The criterion used there was a maximum of |Vfcnfc|, i.e. 
the locus of the steepest drop of n^. This would in fact 
be quite consistent with the present results. However, 
the same caveat as in the present case applies, i.e. this 
criterion will overlook Fermi level crossings of bands with 
low spectral weight [|lO| . 

In our opinion the strange dependence of Vf on elec- 
tron density makes it questionable whether the 'spectral 
weight Fermi surface' in our data is a true constant en- 
ergy contour for a system of 'quasiparticles'. It may be 
possible that at the temperature we are studying a Fermi 
surface in the usual sense no longer exists, and that the 
Hubbard I approximation merely reproduces the spectral 
weight distribution in this case. As our data show, how- 
ever, for that purpose the approximation is considerably 
better than commonly believed. 
Zero temperature studies for the doped t-J and Hubbard 



model are only possible by using exact diagonalization 
]l2[ , in which case the shell-effects due to the small sys- 
tem size require special care fl^jlj] . One crucial point is 
the very different shape of the quasiparticle dispersion at 
zero temperature. Whereas the A band is at least topo- 
logically equivalent to a nearest neighbor hopping disper- 
sion, with minimum at (0,0) and maxiumum at (7r,7r), 
the zero temperature data ||T^ show a second-nearest 
neighbor dispersion with a nearly degenerate band max- 
imum along the antifcrromagnetic zone boundary, and a 
shallow absolute maximum at {tt/2,tt/2). The effect of 
hole doping at zero temperature, however, has a qual- 
itatively very similar effect as in the present case [ [l3[ : 
the chemical potential simply cuts into the quasiparticle 
band for the insulator, which thus is populated by hole- 
like quasiparticles [Q. Again, these 'hole pockets' vio- 
late the Luttinger theorem, indicating again the break- 
down of adiabatic continuity in the low doping regime 
persists also at low temperatures. 
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